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Abstract 

This paper considers the famous Fermi-Pasta-Ulam chain with periodic boundary con- 
ditions and quartic nonhnearities. Due to special resonances and discrete symmetries, 
the Birkhoff normal form of this Hamiltonian system is completely integrable, as was 
shown in |16|. We study how the level sets of the integrals foliate the phase space. 



Our study reveals all the integrable structure in the low energy domain of the chain. 
If the number of particles in the chain is even, then this foliation is singular. The 
method of singular reduction shows that the system has invariant pinched tori and 
monodromy. Monodromy is an obstruction to the existence of global action-angle vari- 
ables. The pinched tori are interpreted as homoclinic and heteroclinic connections 
between travelling waves. Thus we discover a new class of solutions which can be de- 
scribed as direction reversing travelling waves. They remarkably show an interesting 
interaction of normal modes in the absense of energy transfer. These solutions can 
easily be observed numerically. 

1 Introduction 

The Fermi-Pasta-Ulam chain with periodic boundary conditions is a model for point 
masses moving on a circle with nonlinear forces acting between the nearest neighbours. 
The point masses constitute a one dimensional anharmonic lattice. 

This lattice is modelled as follows. Let qj S M be the position of the j-th particle 
(j = 1, . . . ,n) with respect to a certain reference position on the circle. The space of 
positions q = {qi, . . . , g„) of the particles in the chain is M". The space of positions 
and conjugate momenta is the cotangent bundle r*R" of M", the elements of which 
are denoted {q,p) = {qi,--- ,qn,Pi,--- ,Pn)- T'*R" is a symplectic space, endowed 
with the canonical symplectic form dq A dp = ^2^=1 d'lj ^ ^.pj. Any smooth function 
H : r*M" — > M induces the Hamiltonian vector field Xh given by the defining relation 
LXfiidq A dp) = dH. In other words, we have the system of ordinary differential equa- 
tions q, = H-, p, = -f|. 
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The periodic FPU chain with quartic nonhnearities is the special Hamiltonian sys- 
tem on r*R'^ induced by the Hamiltonian 

^= E lp] + Wiq,+i-q,) , (1.1) 

iez/„z 

in which : M — > M is a potential energy function of the form 

W{x) = + \x' . (1.2) 

The system ( |1.1[ ) was introduced by E. Fermi, J. Pasta and S. Ulam Q. These au- 
thors used the chain to model a string of which the elements interact in a nonlinear 
way and they intended to study the ergodic properties of such a string. In their fa- 
mous numerical experiments, they surprisingly observed quasiperiodic behaviour and 
no thermalisation or ergodicity. This puzzling phenomenon finally led to the discov- 
ery of for instance solitons and stimulated the work on the KAM theorem. In his 
recent book |2^, Weissert gives a rather complete overview of the history of the 'FPU 
problem' and its impact on other branches of mathematics and physics. Ford's review 
article [0] contains a nice list of unsolved problems. 

Many different types of explanations were proposed for the observations of Fermi, 



Pasta and Ulam. Nishida |14] and Sanders |18| studied the Birkhoff normal form of 
Hamiltonian ( |1.1D . If this Birkhoff normal form were nondegenerately integrable, then 
the KAM theorem would prove that many low energy solutions of (1.1) are quasiperi- 



odic. These authors were forced to make very restrictive nonresonance assumptions 
though, so they never produced a proof. 



In 1 16 1 the Birkhoff normal form of the periodic FPU Hamiltonian (IT) was finally 
proved to be an integrable approximation of ( |1.1| ). This remarkable result is caused 
by special resonances and discrete symmetries in the Hamiltonian (|l.l| ). The Liouville 
integrable Birkhoff normal form of ( [L.l| ) is the subject of this paper. 

It turns out that the structure of the integrable approximation depends strongly on 
the parity of the number of particles n in the chain. If n is odd, then all the integrals 
of the Birkhoff normal form are quadratic functions of the phase space variables {q,p). 
These quadratics form a set of global action variables. This implies that the foliation 
of the phase space in Liouville tori is trivial and that the equations of motion can be 
solved explicitly. 

The situation is not so simple though, if the number of particles n is even. In this 
case one finds quadratic as well as quartic integrals. They give rise to a singular foli- 
ation of the phase space. We shall find singular fibers that contain relative equilibria 
of focus-focus type. And the singular fibers on which they lie are pinched tori. These 
pinched tori are homoclinic or heteroclinic connections between the focus-focus relative 
equilibria. We explicitly derive under what conditions these connections exist. 

In order to detect the pinched tori and the bifurcations of the relative equilibria, we 
use purely geometric arguments based on invariant theory and the method of singular 
reduction, see This geometric approach reveals all the integrable structure that is 
present in the Birkhoff normal form. It is for example well-known Q, |^] that the 
presence of a pinched torus results in nontrivial monodromy: the fibration of the phase 
space in Liouville tori is not trivial in the sense that the Liouville tori do not form a 
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trivial torus bundle over the set of regular values of the integrals. Instead, we know 
how the Liouville tori are glued together globally, for instance on an energy level set. 
Nontrivial monodromy is an important obstruction to the existence of global action- 
angle variables, see 

We expect that much of the integrable structure of the normal form is still present 
in the low energy domain of the original FPU chain (|1.1|), meaning that the KAM tori 
are 'glued together' in a similar way. 

On the other hand, our analysis yields interesting dynamical information. Previ- 
ously mentioned relative equilibria can be interpreted as waves in the periodic FPU 
chain that travel clockwise and anti-clockwise. Thus it appears that there are ho- 
moclinic and heteroclinic connections between these travelling waves. In the original 
system ( |1 . l]) , we indeed find these direction reversing travelling waves numerically. 
They constitute a large collection of interesting new solutions of the periodic FPU 
chain. These solutions have the remarkable property that they show interesting inter- 
action of normal modes without transferring energy. 



2 Phonons and reduction of the total momen- 
tum 

The frequencies of the linear periodic FPU problem are: 

ujj :=2sin(^) , j = l,... ,n-l . (2.1) 
n 

Note that ujj = oon-j- 

We would like to view the solutions of the Hamiltonian system (|1.1| ) as a superpo- 
sition of waves. Therefore we make the following Fourier transformation: 

Ij = \ )^*: + '^^V )Qn-k + -j^Q^ + ^Qn 2.2) 

\ nujk \ n n ) y/2n " Vn 

Of course, the term with the subscript ^ only appears if n is even. 

The mapping q t-^ q is a linear isomorphism of the position space M". It induces a 
symplectic transformation {q,p) iQ,p) on r*M"'. The Fourier coefficients {q,p) are 
known as 'phonons' or 'quasi-particles'. In phonon-coordinates, the Hamiltonian (|l.l| ) 
reads 

^ n— 1 

^ = + E Y^P' + ^i) + ^4(9-1, • • • ,qn-i) . (2.3) 
J=l 

iJ4 denotes the quartic part of H. An exact expression for i?4 in terms of the ijj can 
be computed by combining ( |1.1| ) and ( |2.2| ). 

The pair of conjugate coordinates {qj,Pj) is called the j-th normal mode. For 1 < 
j < §) we say that the j-th and the (n — j)-th normal mode have wave number j, 
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because both are Fourier coefficients of waves with wave number j, that is waves with 
wave length 1/j times the lenght of the chain. The j-th and (n— j)-th normal mode also 
have the same linear vibrational frequency: ojj = tOn-j- But the j-th and (n— j)-th nor- 
mal mode are out of phase, because one is described by a cosine and the other by a sine. 

Hamiltonian ( |2.3D has an obvious symmetry, namely the flow of Xp^ = ^|-. Clearly, 
this flow defines a free symplectic action on r*R" that leaves the Hamiltonian H in- 
variant. The reduction of the symmetry is standard, cf. Q or by just forgetting 
about the constant ^p^, H trivially reduces to a Hamiltonian function of the variables 
(gi,... ,Pn-i) G T*M""^. Thus we removed the total momentum from 

the equations of motion. 

Because uji, . . . , Wn-i > 0, we conclude using Morse- lemma Q that the origin is a 
stable equilibrium of the Hamiltonian system on T*M""^ induced by ( p. 31 ). It corre- 
sponds to an equidistant configuration of the n particles in the chain. In the sequel we 
will study the flow of the vector field induced by ( p.3| ) on T*R"~^ in a neighbourhood 
of this equilibrium. 



3 The Birkhoff normal form is integrable 



We wish to simplify ( p. 3D by studying its average: 



H{q,p) := hm / F(e*^«2 {q^p)) dt . 

1 —>oc Zl J _'][• 



(3.1) 



The flow of the linear vector field is easily calculated; it is quasiperiodic with 
frequencies loj. 

H is called a Birkhoff normal form for H. Another convenient method to calculate 
normal forms, especially of higher order, is the well-known Lie-series method, described 
for instance in and [l^. In all cases, one knows that there are canonical coordi- 
nates in a neighbourhood of the origin of the phase space r*R"-\ with the following 
property: in these coordinates, H looks like a higher order perturbation of H. 

Thus, under certain nondegeneracy conditions, the low energy solutions of the 
Birkhoff normal form ( |3.lD approximate the low energy solutions of the original system 
(|2.3| ). On the other hand, this is only true for finite n and the domain where the 
approximation is valid probably shrinks when n grows. 

Before presenting H, let us introduce some notation that will be convenient later 
on. For 1 < j < §, define the 'Hopf variables' 



together with 



QjPn—j Qn—jPj 
QjQn~j ~^ PjPn—j 



Pn—ji 



(3.2) 



And if n is even 



^,-2.-2,-2 1-2 \ 
-{qj + Pj + Qn-j + Pn-j) 



Tin := -(on -|- P", 

2 2 2 2 



(3.3) 
(3.4) 



4 



These quantities satisfy the following relations: 



n 



In terms of the Hopf variables, the following expression for H was found in fl^ : 




l<fc<s 



i<fc<f i<fc<f J (3.6) 

In formula (3.6) it is understood that terms with the subscript ^ and j only appear if 
2 respectively 4 divides n. 

As was shown in [^] , due to discrete symmetries and special eigenvalues of Hamil- 
tonian (|1.1|), H remarkably is Liouville integrable. 



Proposition 3.1 // the number of particles n is odd, then the BirkhoJJ normal form 
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(3.6) is Liouville integrable with the quadratic integrals Tij^Uj (1 < j < " ^ 



The integral 7ij has the interpretation of the linear energy of the two modes with wave 
number j, whereas Uj is the angular momentum of these modes. 

Since all these integrals induce a 27r-periodic flow, they constitute a set of global 
action variables. The solutions of the normal form can explicitly be written down. And 
the foliation of the phase space into invariant tori is trivial in the sense that the set of 
Liouville tori is a trivial torus bundle over the set of regular values of the integrals. The 
normal form turns out to be even nondegenerate in the sense of the KAM theorem, 
which proves the abundance of quasiperiodic solutions in the low energy domain of the 



original system (IT). These statements were all proved in [!£]. Thus, this system is 
quite tractable. 

But the situation is not so easy if the number of particles in the FPU chain is even. 



Proposition 3.2 If the number of particles n is even, then the Birkhoff normal form 
(3.6) is Liouville integrable. The integrals are the quadratics 7ij (1 < j < §), Ij ■= 
Uj — u^~j (1 < i < j) and J := ^-^wn (if n is a multiple of A) and the quartics 

■= - ^j^^-j) - - (1 < J < i)- 

The rest of this paper shall be devoted to studying the dynamics and bifurcations of the 
integrable system (3.6) in the case that n is even. We shall encounter nice phenomena 
such as whiskered tori with heteroclinic and homoclinic connections. They shall have 
the interpretation of direction reversing travelling waves and they constitute a class 
of solutions displaying an interesting interaction of Fourier modes with different wave 
numbers in the absence of energy transfer. 



We also study how the level sets of integrals in proposition 3.2 foliate the phase 
space T*M'^~^. This study reveals all the integrable structure that is present in the 
low energy domain of the original FPU chain ( |1.1| ). It turns out that the foliation in 
Liouville tori is not trivial if n > 6 is even. There is monodromy and global action-angle 
variables can not exist. 
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4 Uncoupled equations and regular reduction 



Let us use the short-hand notation J = ^ — j. 

We want to study how the level sets of the integrals of proposition 3.2 foliate the 
phase space T*R"^^. Therefore it is very useful to note that these integrals are uncou- 
pled in the following sense: the integral Tin depends only on the phase space variables 
{qiiL,P!i). The integrals TCil and J depend only on the variables (qn, q3n,,pn,,p3n,). 

2 2 4 4444 

The integrals Ti.j,Ti.j,Ij and ICj depend only on the variables {qj,qj, qn-j,qn-j, PjiPp 

Pn-j ) Pn~j) ■ 

Therefore it is sufficient to know how Tin foliates T*M, how Hn. and J foliate T*M? 
and how Hj^Hj^Ij and ICj foUate r*]R^. The foliation of T*!""^ 
cartesian product of these foliations. 



by all integrals is the 



4.1 The foliation of T^R 

It is clear that Tin = \(q\ + p\) foliates T*M in circles which have their center at the 
origin. In other words, the Fourier mode with the highest wave number waves in 
the chain) and the highest vibrational frequency (wn = 2) has constant energy. 



4.2 The foUation of T* 



The system on T*]R with integrals Tiii and J describes the interaction of the ^-th and 
^-th normal mode. These modes have equal wave number ^ and equal vibrational 
frequency \/2, but are out of phase. 

The foliation of T*M^ can be described in the following standard way. It is clear 



that imHi 



and that for £ 



the inverse image Hn^{hiL) = , . Here 

— 44 , 2hn 



induces an S^-symmetry on 



C T*M is the three-dimensional sphere with radius ^/2hR. The flow of X-)-i„ 

4 ^ 



2h, 



given by 



/ ^5 \ 



qsn 

4 

^ 4 



) 



/ qii COS 621 +piism(j)VL \ 
1^4 ^4 ^4 ^4 \ 

Q3n cos (pn^ + P3n SlU 6 11 

^4 ^4 



cos mn 

^4 ^4 
y p-in COS (f)ll 



^4 ^4 



qsn sm(p2L J 



(4.1) 



The orbits of this circle action have dimension one ii hn > and dimension zero 

4 

otherwise. The reduced phase space can be found by applying the Hopf 

, / 2h n 



map 



2hr 



which is defined by 



^('^ : iq,P) 



(uiL,vn.,wn.) , 

\ 4 4 4 ' ' 



see @. The fibers of J^(^) are exactly the orbits of the S'-'^-action, so S^^ constitutes the 
reduced phase space. Every Hamiltonian function on T*]R^ that commutes with Tiii 
reduces to a Hamiltonian on S? because it is constant on the orbits of the S'^-action. 

H n 
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In particular, J 



1 



2V2n 



. The foliation of S'i in level sets of J is trivial: the level 

Tin ^ 



sets are the circles of constant v^. There are two stable relative equilibria. 




The reduced space S\ „ and the level sets of J 



Reconstructing this picture by the Hopf map, we get the desired foliation of T 



4.3 The foliation of T^R^ 

If n > 6 is even, then for each 1 < j < ^, we have the four commuting integrals 
TLj^TLpXj and K,j on r*M^. They describe the interaction between the two modes with 
wave number j and the two modes with wave number J. Note that these two pairs 
of modes do not interchange energy, since 7ij and 7ij are constants. Still it turns out 
that their interaction is very interesting. The foliation of T*M^ is difficult and contains 
singularities. 

On T*M^, consider the commuting Hamiltonians Tij and Ti^ defined in ( |3.3| ). They 
define a mapping Tijj := {Ti.j,Ti.j). It is easy to see that im Tijj = (M>o)^, whereas 

— modulo permutation of 



/2hi 



for the level sets we know that Tij ^{hj,hj) = S 

coordinates. Here, S| indicates the three-dimensional sphere with radius e. Note that 
e = is not excluded. 

Because the flows of the Hamiltonian vector flelds of TCj and TCj commute and are 
periodic with period 27r, they define a linear symplectic T^-action on x , 



J2h, 

given by 



/ 9j 




\ 




/ Qj COS 


+ pj sin (pj 




\ 


9J 








COS </)J 


+ pj sin (pj 






qn- 


-3 






Qn-j COS 


(pj + pn-jsinc, 


h 




Qn- 


-j 






Qn-j COS 


(pj + pn-j sin 


(Pj 




Pj 








pj COS 


— Qj sin (pj 






Pj 








pj COS (pj 


— p j sin (pj 






Pn 


-5 






Pn- J COS 


(pj - Qn-j sin (, 


h 




\ Pn 


-j 


) 




\ pn-j COS 


(pj - Qn-j sin 


(Pj 


/ 



(4.2) 



The orbits of this T^-action are all tori of dimension #{A; G {j, J} | /ifc 7^ 0}. 



We want to study the reduced phase-space {S 



S^^)/t2. We do this by 



7 



applying a second reduction map 

^(2) : 53 X 53 Six Sl_ defined by 



J^(^) : {q,p) ^ {uj,Vj,Wj,Uj,Vj,Wj) . 

J^^'^^ is the cartesian product of two Hopf mappings. The fibers of J^^^) are exactly 
the orbits of the T^-action, so 5^, x constitutes the reduced phase space. Every 
Hamiltonian function on T*M^ that commutes with TCj and TCj reduces to a Hamiltonian 
on S*^. X Sl because it is constant on the orbits of the T^-action. In particular Ij and 
ICj. The Hamiltonian equations of motion induced by such a Hamiltonian reduce to 
Hamiltonian equations on the orbit space x S^. For a Hamiltonian function H, 
these reduced equations read 

for k = j, J. Thus we obtain a two degree of freedom integrable Hamiltonian system on 
the product of two spheres. We will now study this integrable system, which describes 
the interaction of waves with wave numbers j and j. 





5 Travelling waves 



Consider the integrable system on S*? . x S? 



with Hamiltonian ICj and momentum 



Uj — U7j. The flow of Xj^ induces a symplectic S'^-action on Sf^, x Sf^_ given by 



\ J 



Uj 

Vj COS 2t + Wj sin 2t 
Wj cos 2t — Vj sin 2t 

Uj 

Vj cos 2t — Wj sin 2t 
\ Wj cos 2t + Vj sin 2t j 



(5.1) 



This action has four isolated fixed points, namely the points (it/ij, 0, 0, it/ij, 0, 0). 
But because the Hamiltonian /Cj is invariant under the action (5.1), this implies 
that the derivative of K,j also vanishes at these points. In other words, the points 
(±/ij, 0, 0, ±/ij, 0, 0) constitute the set of joint critical points of Xj and /Cj. 

Critical points of the reduced system on x S\_ are called relative equilibria, 

because in the reconstructed system on T*IR^ - or T*M"~^ if you like - their fibers 
represent invariant sets, see [Q. It follows from (4.2) that in T*R^ the critical fiber 
(.7^'-^))~^(±j/ij, 0, 0, ±j/ij, 0, 0) is the following parametrised torus 



{ {yhj cos ( 



hj cos </>j, y hj sin 0j, 



hj sin (\)j , 



(5.2) 



=FjV/ij sin0J,=FJ 



/lJsin(/>J, ibj 



hj cos ibj 



hj cos (\)j ] 



G } 



It has dimension #{A; G {j,j}\hk 7^ 0}. This torus is invariant under the flow of Xjj so 
one can write the equations induced by H as equations for </) G T^. Using expression 



8 



(3.6) it is not hard to compute that they read: 



d(pj 
~dt 




16 


d(j)j 


dH 


16 



(5.3) 
(5.4) 



Hence the motion in the critical fibers is uniform. The corresponding solutions have a 
clear physical interpretation: with (^]^) one can make the transformation back to the 
original coordinates: 



2/i,- 



cos 



,2Trjk 



+ 



n 



l2h-. .2-Kjk 

COS( (pn) 



In other words, a solution that lies in a critical fiber is a superposition of 

1. a travelling wave = cos(^^^) with amplitude \/^^ and speed approximately 



2. a travelling wave qk = cos(^^^) with amplitude and speed approximately 

The constants hj and hj are supposed to be small, since otherwise the normal form 
approximation has no validity. Thus such a solution is a superposition of two travelling 
waves with wave number j and j. If zizj = ibj then these waves move in the same 
direction. Otherwise one moves clockwise and the other moves anti-clockwise. 

Travelling waves and superposed travelling waves have previously been studied in 
infinite FPU chains [11|. But they can obviously also occur in finite periodic chains. We 
shall study the stability of the travelling wave solutions and homoclinic and heteroclinic 
connections between them. We do this of course in the reduced context, that is we 
consider them as critical points in the reduced phase space 



xSl. 



6 Stability of the relative equilibria 

We want to determine the stability type of the superposed travelling wave solutions in 
the Birkhoff normal form, that is the stability type of the relative equilibria (±/ij, 0, 0, 
±/ij, 0, 0). For this purpose we introduce local coordinates on S*^^ x S^^ near {±hj,0, 0, 
±hj, 0, 0) by projecting {uj,Vj,'Wj,Uj, Vj, Wj) ^ {vj,Wj,Vj, Wj). Note that these are not 
Darboux coordinates. The critical points themselves are all mapped to (0,0,0,0). 



6.1 A Lyapunov function 

One way of proving stability is by pointing out a Lyapunov function. The Hamiltonian 
ICj is the first candidate since it is an a priori constant of motion. But it turns out 
that JCj is not definite at any of the relative equilibria. Luckily, we have another 
constant of motion, namely Ij. We now use that at (itj/ij , 0, 0, ^jhj, 0, 0) one may write 

Uj = ±j^hj-v]-w] = ±j{hj - ^{v] + w]) + ...) and Uj = ±j^h]-v1-w] = 
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which is definite at (0, 0, 0, 0) if and only if itj ^ itj. We conckide that the relative 
equilibria ib(/ij,0,0, — /ij, 0, 0) are stable. 

Near the relative equilibria ib(/ij, 0, 0, /ij, 0, 0) we will try to make linear combina- 
tions of K,j and Xj that are definite. It is easily computed that in local coordinates 

{Vj,Wj,V-j,Wj) 

{32nlCj±2XhjhjIj) = {ujj +Xhj){vj +w'j) + {i^'^ — Xhj){v'^ +w'~)+4ujj(jOj{vjVj—WjWj)+ . . 

modulo constants. Using that \vjVj — WjWj\ < \\{vj,Wj)\\ • \\{vj,Wj)\\, one sees that this 
expression is definite if and only if 



det 



2ujjUJj 



(6.1) 



2 Xh - I ~ ~^ ^^^^ ^^'-^3^3 ~ ^j^j^ ~ ^^j^J ^ *-* 

The preceding inequality has real solutions A if the discriminant 

/I O /I o 

r := ojjhj — lAujjUjjhjhj + ^jhj 
is positive. So if r > then the relative equilibria ±{hj,0, 0, hj, 0, 0) are stable. 

6.2 Linearisation 

Because we still don't know anything about stability if r < 0, an alternative is to 
study the linearisation of the vector field X/Cj at it (/ij, 0, 0, /ij, 0, 0). Again in local 
coordinates, it reads 



1 ° 


-u;]hj 





2u> jUJjhj 


\ 


/ 




\ 







2u> jUJjhj 









Wj 







2ujjUJjhj 





-wf/ij 











\ 2ujjUJjh^ 







/ 


V 


Wj 


/ 



+ h.o.t. 



'h.o.t.' stands of course for 'higher order terms'. One calculates that the characteristic 
polynomial of the above matrix reads 

C(A) = + X'^iiojh] + Lofh] - 8uj]uj]hjh^)/{8nf + 9icojujlh]h])/{8n)^ , 

so the eigenvalues are the numbers 

A 

where 



2uj] hj 



2u;]hj 



/I O AO 

p := 16uJjUJjhjhj — 2ujjhj — 2uJ^h^ , < 

and r as defined previously in ([6.11). Note that C(A) = C(— A), so if A is an eigenvalue 
of ( |6.2| ), then so are —A, A and —A. The reason is of course that our matrix is conjugate 
to an infinitesimally symplectic matrix. 

The next observation is that if r > or g = then p it q^/r G M so the eigenvalues 
are purely real or purely imaginary, dependent on the signs of p ± g-y/r. On the other 
hand, if r < and g 7^ then none of the eigenvalues lies on the real or the imaginary 
axis. A simple analysis now leads to the following results: 
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• If r > then there are four distinct purely imaginary eigenvalues. 

• If r = 0, we find double imaginary eigenvalues. The linearisation matrix is not 
semisimple. 

• If r < then none of the eigenvalues lies on the imaginary axis. 

The set of hj, hj for which r = consists of two half lines in the positive quadrant. 
On these half lines a Hamiltonian Hopf bifurcation [12| occurs: two pairs of imaginary 



eigenvalues come together and split into a quadruple of non-imaginary eigenvalues. 
So the linear stability of the relative equilibria ib(/ij, 0, 0, /ij, 0, 0) changes here from 
neutrally stable to unstable. The linear instability implies of course that the equilibria 
are unstable also in the nonlinear system. This concludes the stability analysis of the 
relative equilibria. 

6.3 Stability results 

We summarize the results of this section in the following 

Corollary 6.1 The relative equilibria ±(/ij, 0, 0, — /ij, 0, 0) are stable. The stability of 
the relative equilibria ib(/ij, 0, 0, /ij, 0, 0) depends on the bifurcation parameter 

r = tOjh'^j - Uu'jujjhjhj + a;|/ij . 

They are stable for r > and unstable for r < 0. 
This is illustrated in the following bifurcation diagram: 




In terms of travelling wave solutions, one may interpret corollary 6.1 as follows: 

The superposition of two travelling waves with wave numbers j and n/2 — j is sta- 
ble if the two waves move in opposite directions. A superposition of waves in equal 
directions can be both stable and unstable. It is stable only if one of the waves is 
relatively small with respect to the other. Otherwise it is unstable. 
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7 Singular reduction 



In the next sections we shall use the method of singular reduction |g] to understand the 
bifurcations of the previous paragraph geometrically, and also to show that if r < 0, 
then there are homoclinic and heteroclinic connections between the unstable relative 
equilibria. 

The flow of Xj^ induces a symplectic S^-action on Sf^, x 5^, given by (5.1). The orbits of 
this 5'^-action are all circles, except for exactly the relative equilibria (it/ij, 0, 0, ±/ij, 0, 0). 
The following quantities are invariant under this action: 

TTj := Uj, pj := Up aj := VjVj — WjWj, Tj := VjWj + v^Wj . 

They satisfy the equation 

The set 

Ph„h, ■■= {(7r,,p„a,-,r,) G MV| + ^ = [h] - T,]){h} - p]), \^,\ < hj, \pj\ < h^} 

therefore constitutes the orbit space Sf^. x S^./s^ Note that we did not yet restrict 
ourselves to orbits of constant Iji this will come later. 

Every function on S^ . x S*^, that commutes with Ij, reduces to a function on Phj,hj 
since it is constant on orbits. In particular, 

= - Pj and ICj := ^(4wjWjcJj - uj^tt] - uj}p]) . 
The reduction map is the map 



jr(^) : {uj , Vj , Wj , Uj, Vj, Wj) ^ {TTj, Pj , aj , Tj j 



which goes from 5"^. x S*^, to Phj,hy The reduction map is a submersion everywhere, 
except of course at the relative equilibria. Unfortunately it is not possible yet to make 
a drawing of Phj,hj since it can not be embedded in M^. But there is an elegant way 
to overcome this problem. 

One sees that both Ij and ICj are independent of Tj, that is these Hamiltonians are 
invariant under the Z2-action generated by the mapping Tj — Tj. This symmetry 
is induced by the discrete symmetries of the original Hamiltonian ( |1.1| ), cf. |16]. The 
orbits of the Z2-action consist of one point if Tj = and otherwise of two points. 
One can reduce the Z2-action by simply forgetting about tj: the reduction map is 
(tTj , Pj , aj ,Tj) ^ {iTj , Pj ,(Tj). The reduced space is the set 

Ph„h-Jz2 = {(^i^Pi^f^j) e MVi < [h] - - Pi), kil < hj, \pj\ < hj} . 

Below, we draw Phj,h-/i2 foi^ hj,hj > : 
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Phj,hj/i2 = ^hj ^ ^h-l S^xi2 shape of a solid pillow. The corners of the pillow 

are cone-like singularities that represent the relative equilibria (it/ij, 0, 0, ±/ij, 0, 0). 

The level sets of Zj = ttj — pj are two dimensional planes. The intersection of such 
a plane with the pillow is a topological disk, a point or empty. Near the singularities 
(vTj, Pj, (Tj) = ±{hj,—hj,0) the disks are very small, indicating that the relative equi- 
libria ±(/ij, 0, 0, —hj, 0, 0) are stable. But near the other two corners of the pillow, the 
singularities it (/ij, /ij, 0), the level set of Ij intersects the pillow in a very large set: 



Tj near hj + /ij Tj near hj — hj 

Let us consider the level set of Ij that passes exactly through the singularity {hj, hj, 0). 
It is the plane Ij = iTj — pj = hj — hj. The intersection of this plane with the pillow is 
the topolocal disk 

{(vTj, cjj) G R^\aj < {hj - TTj){h'j - (vTj -|- hj - hj)"^) , pj = TVj + hj - hj} . 

If hj = hj, then this disk has two singular points. It has one singular point if hj ^ hj. 
The intersection of a level set of ICj with this plane is a parabola. The parabola that 
contains the singularity {hj , hj, 0) is given by the formulas 

aj = a{7Tj) := -^^{^^1 + ^^i^^j + % - hj)^ - ^jh] - u:jh]) , pj = Hj + hj - hj . 

We now make a linear approximation to both this parabola and the singular disk at 
the singular point (/ii,/i2,0). So we calculate the derivative ^ \irj=hj = 2a;^a;- ('^i^i ~'~ 
u}'jhj). On the other hand, the cone-like singularity of the singular disk is approximated 
by 

{|(Tj| < 2\/hjhj{hj — VTj), VTj < hj, Pj = iTj + hj — hj} . 
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So the tangent hne to the parabola points into the cone exactly if —2^yhjhj < 2Z7^('^j 
iOjhj) < 2^yhjhj, that is if r < 0. In this case, the critical point (-Kj, pj, aj) = (hj, /ij, 0) 
is clearly unstable, which agrees with our previous analysis. The tangent to the 
parabola does not point into the cone if r > 0. The following two pictures repre- 
sent the two possibilities: 



r > r < 

In the case that r > 0, we now see that ICj, restricted to [Phj,h-Ji2) ^1J^{hj — %), 
is extremal at {hj,hj,0). In other words, (/ij, 0, 0, /ij, 0, 0) is a stable relative equi- 
librium inside the singular fiber (5"^^. x S"^,) riT~^{hj — hj). A theorem of Montaldi 

1 13] then states that {hj,0,0,hj,0,0) is stable in S"^, x Sf^_. This is an entirely geo- 
metric way to see the Hamiltonian Hopf bifurcation and the stability change of the 
relative equilibrium (/ij, 0, 0, /ij, 0, 0). The same argument holds for the equilibrium 
(-/ij,0,0,-/ij,0,0). 



8 Pinched tori and monodromy 

We shall argue that the reduced system on S^. x S^, induced by ICj has nontrivial 
monodromy and that there are homoclinic and heteroclinic connections between the 
relative equilibria ±(/ij, 0, 0, /ij, 0, 0) if r < 0. 

This can all be shown by simply drawing a picture of {Phj,h-Ji2) ^1J^{hj — hj) 
and its foliation into the level sets of K,j. Recall that {Phj,hj/i.2) ^1J^{hj — hj) is a 
topological disk that contains one singular point if hj ^ hj and two singular points if 
hj = hj. Each point inside the disk represents two three-dimensional tori in r*M^. The 
regular points on the boundary of the disk each represent one three dimensional torus. 
The singular points, which also lie on the boundary of the disk, represent a 'singular' 
two-dimensional torus. This singular torus has the interpretation of a superposition of 
travelling waves in the FPU chain. 

Let us first consider the case that hj = hj. In that case the reduced space is 
bounded by parabolas 

nTri(o) = {{-Kj^pj^aj) eW^\\aj\< h] - , pj = vr^} 

It has two singular points which lie on the same level set of ICj . This level set is also a 
parabola. Furthermore, note that r = {(jjj + coj — 14u;|a;|)/i| = 16(1 — uj'juj'j)h'j, which 
is negative if and only if ujjCOj > 1 if and only if n/12 < j < n/4. But then the reduced 
space simply looks like this: 
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We observe that there is a heterochnic connection between the two singular points. 
In the reconstructed system on S"^. x S*^, this corresponds to a doubly pinched torus: 
two focus-focus singular points of which the stable and unstable manifolds coincide. 
In the original phase space T*R^, this doubly pinched torus is again reconstructed as 
a heterochnic connection between two-dimensional tori. They are connected by their 
'whiskers' which have dimension four and are both diffeomorphic to M x T^. 

The heterochnic connection has the following interpretation. If one starts with a 
motion that is nearly the superposition of two travelling waves with equal energy and 
in equal direction, then after a certain 'incubation period' both travelling waves will 
come to a halt and turn around until the motion looks very much like the superposition 
of two travelling waves the directions and energies of which are again equal, although 
the direction is opposite to the direction in the beginning. This process continues and 
the superposed waves keep changing direction. 

In the case that hj ^ hj and r < 0, the singular level set generically looks like this: 



We see that there is a homoclinic connection that connects the singular point to itself. 
In the reconstructed system on 5^^. x 5*^, this corresponds to a singly pinched torus: 
the stable and unstable manifold of the focus-focus singular point coincide. In T*R^ 
the pinched torus is a homoclinic connection of a two-dimensional torus to itself. The 
two-dimensional torus again represents the superposition of travelling waves in the 
same direction with wave number j and J, which now do not have equal energy. 

This completely describes the interaction between travelling waves with wave num- 
ber j and travelling waves with wave number n/2 — j. The travelling waves do not 
interchange any energy. Still, their influence on one another is such that each of the 
travelling waves drastically changes its momentum and thus its direction. 

8.1 Monodromy 

It was shown in [^] that the presence of a pinched torus implies monodromy: the 
regular Liouville tori in S?, x S? do not constitute a trivial torus bundle. Instead, 
there is monodromy and the monodromy map is known. When we reconstruct, we 
see that the regular tori in T*M^ and T*W^~^ can not form a trivial bundle either. 
Nontrivial monodromy is an important obstruction to the existence of global action- 
angle variables. Recall that the pinched torus is present if n > 6 is even. We conclude 
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that if n > 6 is even, then the integrable normal form (3.6) does not admit global 
action-angle variables. 

We have now described how the level sets of the integrals of the Birkhoff normal form 
(3.6) globally foliate the phase space. The same foliation describes all the integrable 
structure that is present in the low energy domain of the original FPU Hamiltonian 
( |1.1| ). In a certain sense, the KAM tori in the low energy domain of the phase space 
will also form a nontrivial bundle. 



9 Numerical solutions for 16 and 32 particles 

I tried to detect the direction reversing waves numerically in the original periodic 
FPU system ( |1.1[ ). This system is not integrable, but in the low energy domain it is 
approximated by the integrable normal form. Therefore we expect that the low energy 
solutions of FPU Hamiltonian ( |1.1| ) behave as predicted by the Birkhoff normal form. 

Let us study the periodic FPU chain with n = 16 particles, a number chosen by 
Fermi, Pasta and Ulam themselves. We shall investigate solutions that start out as a 
superposition of two travelling waves with wave numbers j = 3 and j = ^ — 3 = 5. 
Note that = | < j = 3 < ^ = 4, so if we give both waves equal energy and equal 
direction, then we expect the solution to reverse its direction. 

Let us take the initial conditions q{0) = (0, 0, .25, 0, .25, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0) and 
p(0) = (0,0,0,0,0,0,0,0,0,0, .25,0, .25,0,0), such that ^3(0) = ^5(0) = ^3(0) = 
1(5(0) = .0625. The total energy is H = 0.1742. The angular momenta and U5 
measure the direction of the waves. Their values in the course of time are plotted here: 




Superposed travelling waves with wave numbers 3 and 5 in the FPU chain with 16 particles 
interact as predicted by the normal form: the angular momenta exhibit a relaxation 
oscillation between their maximal and minimal values. 

Note that this line is double: the angular momenta remain equal all the time, just 
as was predicted by the normal form. Furthermore, they vary between their maximal 
and minimal values at the given energy. The heteroclinic connection of the normal 
form is clearly visible here: the angular momenta seem to stick to their maximal and 
minimal value for quite some time, before moving off again. It is not very likely though 
that this heteroclinic connection really exists in the original system. Generically, coin- 
ciding stable and unstable manifolds start having transverse intersections under small 
perturbations, resulting in chaotic orbits. On the other hand, the angle under which 
they intersect is extremely small, such that this chaos is hardly visible. Hence our 
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observations do not look chaotic at all. 

To compare, I also investigated how waves with wave numbers 2 and 5 inter- 
act. Note that 2/5. So I took q{0) = (0, .25, 0, 0, .25, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0) and 
p{0) = (0,0,0,0,0,0,0,0,0,0, .25,0,0, .25,0), such that 7^2(0) = 7^5(0) = U2I0) = 
^5(0) = -0625. The total energy is H = 0.1528. The values of U2 and are depicted 
here on the same time scale: 

o.oe| 1 1 1 1 1 1 1 1 1 1 

0.07 - 
0.06 - 



0.02 
0.01 



As predicted, waves with wave numbers 2 and 5 do not interact. 

As predicted by the normal form, waves with wave numbers 2 and 5 do not interact at 
aU. 

1 also investigated how the waves with wave numbers 3 and 5 interact at a higher 
energy level. So we start with q{0) = (0, 0, .5, 0, .5, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0) and p(0) = 
(0,0,0,0,0,0,0,0,0,0, .5,0, .5,0,0), such that 7^2(0) = 7^5(0) = ^2(0) = 115(0) = .25 
and the total energy is H = 0.7048. The result: 




100 200 300 400 500 600 700 800 900 1000 



Also at a higher energy level, the normal form is a good approximation. 

Note that even at this high energy level the normal form still constitutes a very good 
approximation! This is a remarkable observation, but similar phenomena have been 
observed in other systems. 

Finally, 1 integrated the chain with 32 particles with initial conditions qj (0) = pj (0) = 
for aU j except 57(0) = 99(0) = ^23(0) = ^25(0) = .25 such that HviO) = HgiO) = 
U7{0) = ■U9(0) = .0625 and H = 0.1762. The results are still very nice: 
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For the chain with 32 particles, the normal form works perfectly at this energy level. 

But at a higher energy the normal form is no longer a good approximation when 
n = 32. If one chooses the initial conditions gj(0) = Pj(0) = for all j except 
^7(0) = ^9(0) = P23(0) = P25(0) = .5, then ^7(0) = 7^9(0) = n7(0) = ^^9(0) = .25 and 
H = 0.7091. The angular momenta interact as follows: 




At high energy, the normal form does not make much sense anymore for n = 32. 

Thus, the domain of validity of the normal form shrinks if n grows. This is to be 
expexted, because when n is very large, near-resonances start spoiling the validity of 
the normal form. Still, the normal form gives us some surprisingly good predictions, 
even at energies where one would not expect this anymore. 



10 Conclusions and discussion 

We considered the Birkhoff normal form of the periodic Fermi-Pasta-Ulam chain with 
periodic boundary conditions and quartic nonlinearities. In the domain of the phase 
space where the solutions have low energy, this Birkhoff normal form constitutes an 
approximation of the original Fermi-Pasta-Ulam chain. And because of symmetries 



and nonresonance properties, it is Liouville integrable, see |16]. We studied how the 
level sets of the integrals foliate the phase space, because one expects that a lot of the 
integrable structure of the normal form is still present in the low energy domain of the 
original FPU chain. 

The most interesting phenomena happen when the number of particles n is even. In 
that case the integrals are quadratic and quartic functions of the phase space variables. 
Using the purely geometric methods of regular and singular reduction we showed 
that the foliation has singular elements, pinched tori. 



It is well-known P], |19] that pinched tori imply monodromy: the Liouville tori do 



not form a trivial torus bundle over the set of regular values of the integrals. This 
is important information about how the Liouville tori are glued together globally, for 
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instance on an energy level set. Among others, monodromy is an important topologi- 
cal obstruction to the existence of global action-angle variables. Because the Birkhoff 
normal form approximates the FPU chain especially well in the low energy domain, we 
expect that the KAM tori on the low energy level sets are 'glued together' similarly. 

At the same time, our study unravels interesting dynamical information. We are 
able to determine how waves with different wave numbers interact. It turns out that 
waves with wave number j can only interact with waves of which the wave number is 
n/2 — j. And even though these waves do not interchange any energy, their interac- 
tion is far from trivial. The pinched tori that were mentioned before are homoclinic 
and heteroclinic connections between solutions which are a superposition of travelling 
waves with these wave numbers. Thus it can happen that these superposed travelling 
waves change their direction. 

In the original FPU chain (1.1), we indeed find these direction reversing travelling 
waves numerically in the chains with 16 and 32 particles. They form a class of inter- 
esting new solutions of the periodic FPU chain. 

Surprisingly, one can even observe direction reversing travelling waves at rather 
high energy levels, where the Birkhoff normal form is a very questionnable approxi- 
mation. It would be interesting to study how robust the Liouville tori near a pinched 
torus are under Hamiltonian perturbations. Maybe the KAM theorem can produce 
extra strong conclusions for systems with monodromy, which enables us to understand 
the validity of the normal form at high energy. 

Finally, the reader should be aware of other wave reversal phenomena that have 
been observed in the literature. I especially refer to |^ which studies the 'boomeron', 
a soliton that comes back. Thus, we have found yet another interesting link between 
the Fermi-Pasta-Ulam chain and integrable wave equations. 
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